/**
  CNOK project, Anyang Normal University, IMP-CAS
  \brief Calculate inclusive Rs=sigma_exp/sigma_th for -N reactions.
  \author SUN Yazhou, asia.rabbit@163.com
  \since 2022/05/23
  \date Last modified: 2022/05/25 by SUN Yazhou
  \copyright 2020-2023 SUN Yazhou
  \copyright CNOK project, Anyang Normal University, IMP-CAS
*/

#include <vector>
#include "TANucleus.h"

using std::vector;

static const int FNNPAR = 1; // 0: Horiuchi; 1: LenziRay
static const bool PAULI = false; // 0: off; 1: on
// static const double EK = 2100.; // incident energy: MeV/nucleon

int main(int argc, char const *argv[]){
  vector<TANucleus *> ov;
  TANucleus *t;

  // C14 //
  // (A, sp, sn, cs, dcs)
  ov.push_back(t = new TANucleus(14., 8.1764, 20.8310, 80.19, 7.39));
  t->SetBaseDir("config/C/C14");
  t->SetOrbit({"0p1/2_1/2-", "0p3/2_3/2-", "1s1/2_1/2+", "0d5/2_5/2+"});
  t->SetEx({0., 3.6845, 3.0894, 3.8538});
  t->SetC2S({1.58, 2.301, 0.021, 0.118});

  for(auto &t : ov){
    // t->SetPauli(PAULI);
    // t->SetEk(EK);
    // t->SetFNNParOpt(FNNPAR);

    // if(t->GetA() == 16.)
      t->Rs();
  } // end for

  return 0;
} // end the main function
